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Abstract. We consider the Neumann Laplacian with constant magnetic field 
on a regular domain. Let B be the strength of the magnetic field, and let Ai (B) 
be the first eigenvalue of the magnetic Neumann Laplacian on the domain. It 
is proved that B i—* \i(B) is monotone increasing for large B. Combined with 
the results of FoHc2 , this implies that all the 'third' critical fields for strongly 
Type II superconductors coincide. 



1. Introduction and main result 

Let ft C K 2 be a bounded, simply connected domain with regular boundary. We 
keep this assumption in the entire paper. 

Let F(x) = (FijFs) — (—X2/2,x%/2) — a standard choice for a vector potential 
generating a unit magnetic field: curlF = 1. We consider Tt(B), the self-adjoint 
operator associated with the closed, symmetric quadratic form, 

W 1 ' 2 (CI)bu^Q b (u) = [ |(-iV + B¥)u\ 2 dx. 

Jn 

We will use the notation p_\ = (— iV + A). Then, more explicitly, H-(B) is the 
differential operator p 2 BT with domain {u € W 2 ' 2 (Vl) v ■ pbfu\oq — 0}, where v is 
the unit interior normal to dtt. 

We choose and fix a smooth parametrization 7 : ^rS 1 dQ of the boundary. 
We may assume that |7'(s)| = 1 for all s. We will further parametrize -^r-S 1 by 
[— |<3Q|/2, |9f2|/2] with periodicity being tacitly understood. 

For a point p = 7(5) £ dfl we define k(p) — also denoted by k(s) — to be the 
curvature of the boundary at the point 7(5), i.e. 

</'(*) = k(s)u(s), 

where v{s) is the interior normal (to the boundary) vector at the point 7(3). The 
maximum of k will play an important role, we define therefore, fc max := max s {fc(s)}. 

Define \\{B) = inf Spec H(B) to be the lowest eigenvalue of Tt(B). The dia- 
magnetic inequality tells us that 

Ai(B) > A!(0), 
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for all B > 0. 

One may ask whether the more general inequality 

0<B 1 <B 2 => Ai(Si) < Ai(Sa), 

which one can consider as a strong form of diamagnetism, holds (see [Erdlj . |Erd2j 
and |LoThaj ). 

In this paper we prove that strong diamagnetism holds for sufficiently large B. 

Theorem 1.1. 

The one sided derivatives, 

a; +i B) = hm + a; {b) = hm 

exisi for all B > and A' x satisfies 

liminf A; ,(£) > 0. (1.1) 

B-^oc ' 

Furthermore, there exists a universal constant Oo > smc/i i/ia£ i/ f2 is not a disc, 
then the limit actually exists and satisfies, 

hm X' h _(B) = hm X' h+ (B) = O . (1.2) 

If is a disc, then 

limsupA^ + (B) > 8q, 
< liminf A; + (B) < 9 . 

B^oo 

In particular, in any case, there exists Bq > suc/i that Xi(B) is strictly increasing 
on [Bq, oo). 

Results similar to (jl.ll) have been proved recently in [FoHe2j under extra assump- 
tions. First of all (in |FoHel| ) a complete asymptotics of Ai(-B) was derived for ft 
satisfying a certain 'generic' assumption, i.e. that the boundary curvature only has 
a finite number of maxima, all being non-degenerate. This complete asymptotics 
was then used to obtain (jl.lj) . The most prominent domain excluded in this ap- 
proach is the disc — where the curvature is constant. However, |FoH e2 includes a 
special analysis of the disc proving that Theorem 11.11 remains true in that case. 

What remained was the study of all the other 'non-generic' cases. Also it seemed 
desirable to be able to establish Theorem 11.11 without using the existence of a 
complete asymptotic expansion, since such expansions are difficult to obtain and 
their structure depends heavily on the different kinds of maxima of the boundary 
curvature. In this paper we realize such a strategy. It turns out that for all domains, 
except the disc, one can modify the approach from FoHe2j to obtain (jl.lj) with only 
very limited knowledge on the asymptotic behavior of \±(B). For the disc one can 
use the special symmetry (separation of variables) of the domain to conclude. 

Thus the structure of the proof of Theorem 1 1.1 1 is as follows. The statements for 
the disc follow from the analysis in |FoHe2] which will not be repeated. Thus we 
only consider the case where SI is not a disc. If ft is not a disc then there exists 
a part of the boundary where the ground state will be very small. Thus one can 
choose a gauge such that |A^| <C 1 (for large B and in the L 2 -sense), where A is 
the vector field F in the new gauge. This new input to the proof in [FoHe2j allows 
us to differentiate the leading order asymptotics for \i(B). 
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Notice that if f2 is not a disc, then it satisfies the following assumption : 
Assumption 1.2. 

If we denote by II the set of maxima for the curvature, i.e. 

II = {p G dfl | k(p) = fc max }, 

then 

n ^ an . 

Finally, we will prove in Sectionl3l (Theorem 13. 3j) that all the natural definitions 
of the third critical field appearing in the theory of superconductivity coincide with- 
out any other geometric assumption than regularity and simply connectedness. 



2. The analysis of the diamagnetism 

Two universal constants 0o,Ci will play an important role in this paper, as in 
any investigation of the magnetic Neumann Laplacian. For detailed information 
about these constants, one can refer to 1 le.Moj . For the second constant Ci, we 
only use the fact that it is strictly positive. The first, Qo can be defined as the 
ground state energy of the magnetic Neumann Laplacian with unit magnetic field 
in the case of the half-plane, 

6 :=Ai(B = l), for fl = R 2 + . 

The numerical value of Go can be calculated with precision (Oo ~ 0.59), however 
for our purposes the following (easily established) rigorous bounds 

< 9 < 1, 

suffice. 

We recall the following general, leading order asymptotics of Ai(-B) proved in 
[HeMoj . 

Proposition 2.1. 

As B — > +oo, then 

Ai(-B) = 60-B + o(B) . (2.1) 

If a state u is localized away from the boundary, i.e. u <E C^°(fl), we have the 
following standard inequality 

(u, H(B)u) > B\\uf L2m , 

where, for v, w in L 2 (Q), (v , w) denotes the L 2 scalar product of v and w. 
Using that Oo < 1 it is therefore a standard consequence of (|2.ip (for the proof see 
[HeMo] ) that ground states are exponentially localized near the boundary. 

Lemma 2.2 (Normal Agmon estimates). 

There exists a, M, C > such that if B > 1 and ipi{ • B) is a ground state ofH(B) 
then 

f e^dist^an) { ] ^ {x , B) |2 + 1 |psp ^ { .- B )\ 2 }dx 
isi B 

< G I \ip!{x; B)\ 2 dx. (2.2) 

-'{v / Bdist(a;,On)<M} 
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In particular, for all N > 0, 

dist(z, 90)^1^1 (x; B)| 2 = ©(B-^ 2 ). (2.3) 



^From [HcMo, Proposition 10.5] we also get the following (stronger than (|2.ip ) 
result, 

Proposition 2.3. 

Let Qq, C\ be the usual universal constants and define, for C > 



U B (x) 



B, d\st{x,d^)>2B- 1 / & , 
O B - dk(s)VB - CB 1 ' 3 , dist(x, dfl) < 2B~^ 6 . 

Then, if B > 1 and C is sufficiently big, we have for all tp £ W 2 ' 2 (fl), 
Oj,H(B)tP) > f U B (x)\iP(x)\ 2 dx. 



Proposition I2.3I and a corresponding improved upper bound (also proved in 
[HeMop . 

Ai(S) = Q a B-C l k m33! VB + o{VB), (2.4) 

imply, by suitable Agmon estimates, that ground states have to be localized near 
the set n. We actually only need the following very weak version of this localization. 

Lemma 2.4. 

Let eo > 0. Then, for all N > 0, there exists C > such that if ■ ; B) is a 
ground state for Ti(B), then 



[ \in(x; B)\ 2 dx < C B- A 

J {dist(x,n)>e a } 



We now introduce adapted coordinates near the boundary. Define, for to > 
$ : x (0, t ) -» $(s, t) = 7 (s) + tu{s). 

Z7T 

For to sufficiently small we have that dist($(s, t), dfl) = t and that $ is a diffeomor- 
phism with image {ir £ f2 1 dist(a;, dQ) < to}. Furthermore, the Jacobian satisfies 
\D$\ = 1 - tk(s). 

Lemma 2.5. 

Let us define for e < min(to/2, |cT2|/2) and sq G dQ 

O(e,s ) :={a; = $(s,f)|t<e,|s-so|>e}. 
Then there exists (f> 6 C°° (11) suc/i that A = F + V0 satisfies 

|A(x)| < Cdist(a;,9fi), 

/or a; g Sl(e, so). 
Proof. 

Let A = (Ai,A 2 ) be the magnetic 1-form pulled back to (s,t) coordinates, 

F x dx + F 2 dy = Axds + A 2 dt. 
Taking the exterior derivative, and using dx A dy = \D&\ds A dt, we find 
curl s , t A = d s A 2 - d t Ai = (1 - tk(s)). 
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Since {(s,t) \ t < e, \s — sq\ > e} is simply connected there exists a function <f> € 
C DO ($- 1 (n(e,s ))) such that 

A + V M 0= (i - i 2 A:( 5 )/2,0). 

Let X eC°°(n), 

X = l on {x 1 i < e, |s - s | > e}, 

X = on {2; | dist(x, <9ft) > 2e or |s - s | < e/2}, 

and define 0(x) = ^($ _1 (a;))x(a;)- Then solves the problem. □ 

Proof of Theorem ll.il 

Let % G C°°(Ti) be such that F := F + satisfies F • v = on dfl. The existence 
of such a. (j> is easy to prove. Define H(B) to be the self-adjoint operator associated 
to the closed quadratic form 

W 1)2 (fi) 3 «h / | -iVu + BFupcte. 

Then Tt(B) and TL{B) are unitarily equivalent and so they have the same spectrum. 
Furthermore, the domain of H(B) is 

V(H(B)) = {ue W 2a {9) : v • Vu | fln = 0}, 

in particular, T>(7i(B)) is independent of £?. Applying analytic perturbation theory 
to Ti(B) we get the existence of A^ + (B), \[ - (B). 

We recall that Theorem 11.11 was proved already in FoHeT] in the case of the 
disk, so it remains to consider the case where is not the disc. Thus f2 satis- 
fies Assumption 11.21 Therefore, there exist sq € [— |9fi|/2, |9fi|/2] and < eo < 
min(i /2, |dft|/4) such that 

[so-2eo,so + 2e o ]nn = 0. 
Let A be the vector potential defined in Lemma l2.51 Qb the quadratic form 

W l > 2 {n) 3 u h-> Q B (u) = [ | -iVu + BAu\ 2 dx, 

Jn 

and T-t(B) be the associated operator. Then Tt(B) and H(B) are unitarily equiva- 
lent: H(B) = e lB ^Ti.(B)e^ lB ^ , for some <fi independent of B. By analytic pertur- 
bation theory applied to Ti-(B) there exists an analytic branch of eigenfunctions, 

H((3)^ + (-;f3) = Ai •;/?), 

for [3 S [B,B + e), some e > 0, with HV^+C^H = 1. 

With ipt( ' i ft) := e^Vi + (-; (3) being the corresponding eigenfunctions of TC(f3), 
we get 

= (A>f( • ; B) , p BjL i>t( ■;B)} + {p B ^t{ -,B), M+( ■ ; B)) 

+ 2K{Q B (i; ) ^+ (£))}, (2.5) 
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where v = jp^t (P)\ g- B - The last term in (|2.5p vanishes because ipf is a normal- 
ized eigenfunction of TL, and therefore, 

A' 1)+ (B) = (A>+ ( • ; B) , ( • ; B)) + ( Pb ^j+ ( • ; B) , A>+ ( ■ ; B)) . (2.6) 

We now obtain for any (3 > 0, 

a; + = QB +g (^(-^)) -QB(^(-^)) _ ^ / |i|2 ,^ (a;;S)|2 dx 

P JO 

P Jo 
By Lemma |2~51 we can estimate 



f |A| 2 \i>f{x;B)\ 2 dx<C ( &ist(x,dtt) 2 \^+(x;B)\ 2 dx 
Jn Jn 

+ [ \iPt(x;B)\ 2 dx. (2.8) 

Jo\0(£ ,so) 

Combining Lemmas 12.21 andl 2. 41 we therefore find the existence of a constant C > 
such that : 

/ \A\ 2 \ifj+(x;B)\ 2 dx < CB- 1 . (2.9) 
Jo 



We now choose (3 — r\B, where r\ > is arbitrary. By the weak asymptotics (|2.ip 
for Ai(-B), we therefore find : 

limiuf X' 1+ (B) > 6 -ryC. (2.10) 

Since 77 was arbitrary this implies 

liminf A' 1+ (B) > 9 . (2.11) 

Applying the same argument to the derivative from the left, X[ -(B), we get (the 
inequality gets turned since /3 < 0) 

limsupAi < 9 . (2.12) 
Since, by perturbation theory, A^ + (B) < X[ _(B) for all B, we get (|1.2p . □ 



3. Application to superconductivity 

As appeared from the works of Bernoff-Sternberg [BeStj . Del Pino- Felmer- Stern- 
berg [dPiFeStj . Lu-Pan [LuPall lLuPa2l ILuPa3j . and Helffer-Pan [HePaj . the deter- 
mination of the lowest eigenvalues of the magnetic Schrodinger operator is crucial 
for a detailed description of the nucleation of superconductivity (on the boundary) 
for superconductors of Type II and for accurate estimates of the critical field Hc 3 ■ 
In this section we will clarify the relation between the different definitions of criti- 
cal fields considered in the mathematical or physical literature and all supposed to 
describe the same quantity. This is a continuation and an improvement of |FoHe2j : 
we will be indeed able to eliminate all the geometric assumptions of that paper. 
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We recall that the Ginzburg-Landau functional is given by 

£[1>, A] = S K , H [^, A] = J {\ Pk ha^\ 2 - k 2 H\ 2 + y IV'I 4 

+ K 2 iJ 2 |curlA- l| 2 |da; , (3.1) 

with (V>,A) e W^ 2 {n ; C) x W 1,2 (Q,; M 2 ). 

We fix the choice of gauge by imposing that 

div A = in O , A • v = on dn . (3.2) 

We recall that the domains Q are assumed to be smooth, bounded and simply- 
connected and refer the reader to [BonJ , BonDa] and |BonFo| for the analysis of the 
case with corners. 

By variation around a minimum for £ Kj h we find that minimizers (ip, A) satisfy 
the Ginzburg-Landau equations, 

J*HAtf = « a (l-W> 1 in fl- (33a) 

curl^-^W^-^V^-I^A J 

curlA-l = } ° n 90 < (3 ' 3b) 

with 

curl 2 A = (9^ (curl A) , — <9 xi (curl A)) . 

It is known that, for given values of the parameters k, H, the functional £ has 
(possibly non-unique) minimizers. However, after some analysis of the functional, 
one finds (see (GiPhj for details) that, for any n > 0, there exists H (k) such that if 
H > H(k) then (0, Fn) is the only minimizer of £ Kj h (up to change of gauge). 
Here we choose Fn as the unique solution in £1 of cutIFq = 1 satisfying (|3.2p . 
Following Lu and Pan |LuPalj . one can therefore first define 

H_c 3 ( K ) = inf{# > : (0, Fn) is a minimizer of £ k ,h} ■ (3.4) 

In the physical interpretation of a minimizer A), IVK 2 -)! i s a measure of the 
superconducting properties of the material near the point x. Therefore, H_c 3 ( K ) ^ s 
the value of the external magnetic field, H , at which the material loses its super- 
conductivity completely. 

Actually, as already used implicitly in [LuPal and more explicitly in |FoHe2j . 
we should also introduce an upper critical field, H_q (k) < Hc 3 {k), by 

Hc 3 {k) = infjiJ > : for all H' > H, (0, Fn) is the only minimizer of £ K h'} ■ 

' (3-5) 

The physical idea of a sharp transition from the superconducting to the normal 
state, requires the different definitions of the critical field to coincide. 

Most works analyzing H_c 3 relate (more or less implicitly) these global critical 
fields to local ones given purely in terms of spectral data of the magnetic Schrodinger 
operator H(B), i.e. in terms of a linear problem. The local fields are defined as 
follows. 

HcI(k) = inf{ff > : for all H ' > H, Xi(kH') > k 2 } , 

H}§*(k) = inf{H > : X^kH) > k 2 } . (3.6) 
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The difference between Hq (k) and H^(k) — and also between Hc 3 (k) and 
—C 3 ( K ) — can be retraced to the general non-existence of an inverse to the function 
B i ► Ai(S), i.e. to lack of strict monotonicity of Ai. In the previous section, 
we have solved this monotonicity question and we now explain, following mainly 
[FoHe2| . how this permits to close the discussion about this 'third' critical field in 
the high k regime. 

The next theorem, which is proved in [FoHc2 , is typical of Type II materials, in 
the sense that it is only valid for large values of k. 

Theorem 3.1. 

There exists a constant kq > such that, for k > kq, we have 

K C3 («) = H% c 3 (k) , H C3 («) = H l °l( K ) . (3.7) 

On the other hand, we have from Theorem ll.il : 
Proposition 3.2. 

There exists kq such that, if k > kq, then the equation for H : 

\ 1 (kH) = k 2 , (3.8) 

has a unique solution H(k). 

In other words, for large k, the upper and lower local fields, defined in (13. 6| . 
coincide. We define, for k > ko, the local critical field Hq^(k) to be the solution 
given by Proposition 13. 21 i.e. 

Ai(kH£, c (k)) = k 2 . (3.9) 

Using Proposition 13. 21 we can identify the lower and upper local fields and there- 
fore find the following result. 

Theorem 3.3. 

Suppose Q is smooth, bounded and simply connected. There exists no > such that, 
when k > no, then 

H x S:(k) = H C3 (k) =HcM ■ (3.10) 

Remark 3.4. 

This result was established in |FoHe2| under the additional assumption that fl was 
either a disk or a domain whose boundary has only a finite number of points of 
maximal curvature (with in addition some non degeneracy condition). 
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